Importance of on-site interaction in graphene. 
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We use the Gutzwiller method to investigate the importance of the on-site Coulomb interaction 
in graphene. We apply it to Hubbard Hamiltonian to study the renormalization of the kinetic 
energy in graphene due to the on-site Coulomb interaction. We find that a reasonable strength of 
the interaction has a very weak effect in reducing the kinetic energy. Hence we predict that the 
Brinkmann-Rice metal-insulator transition in graphene is not possible. The effect is understood in 
terms of the high kinetic energy in graphene. 

PACS numbers: 



Graphene is a two dimensional (2D) crystal of car- 
bon atoms. It is an allotropy form of Carbon in which 
the atoms form a honeycomb lattice. The experimen- 
tal success of creating the graphene [3, Q has renewed 
the interest of studying the electronic properties of the 
two dimensional electron system. It is because the 2D 
crystals were predicted to be thermodynamically unsta- 
ble at finite temperature. [H, 0] Graphene is stabilized 
by strong covalent bonding between carbon atoms and 
a small crumpling in the third dimension which increases 
the elastic energy. It was known theoretically that the 
low energy excitations of graphene behave as massless 
Dirac fermion Q ; which leads to different electronic prop- 
erties in this system compared to that of the conventional 
two dimensional electron system. In fact, the observation 
of the half integer quantum Hall effect in graphene is one 
of the examples. Q 

In graphene the nearest neighbor hopping, which sets 
a scale to the kinetic energy, is very high. In such case 
we need to understand how important the role of the 
Coulomb interaction will be in determining the electronic 
properties. Our previous study shows that the effect of 
the long range Coulomb interaction is not important in 
single layer graphene. We have shown that this system 
can not have charge inhomogeneous states (@,0|)- It has 
also been shown that this system can not undergo a fer- 
romagnetic transition(Q). Even if there is a consensual 
agrement about the fact that the Coulomb interaction 
does not have a strong effect on the electronic proper- 
ties of graphene, we do here an explicit calculation to 
benchmark this statement. We study the importance of 
the on-site Coulomb interaction using Hubbard model in 
Gutzwiller approach. [9lll0] In Gutzwiller frame work the 
Hubbard model leads to the possibility of a metal insu- 
lator transition at half filling, which is also known as the 
Brinkmann-Rice transition. The idea is to look at the 
effect of the on-site Coulomb repulsion on the renormal- 
ization factor of the kinetic energy which is a direct mea- 
sure of the correlation effect. In graphene we find that 
the renormalization in the kinetic energy is very small. 
Hence we infer that the system can not be even close to 



the Brinkmann-Rice Metal-Insulator transition. 

First, we introduce the Gutzwiller method briefly for 
which we closely follow Ref. [ll[ ■ Among numerous the- 
oretical approaches, the Gutzwiller method provides a 
transparent physical interpretation in terms of the atomic 
configurations of a given site. Originally, it was applied 
to the one-band Hubbard model Hamiltonian. [12| In 
graphene also the electronic band close to the Fermi 
level is formed by the hybridization of neighboring p z 
orbital. So introducing single Slater Koster [13[ parame- 
ter ppir ~ 3.0eV it is straight forward to use the method 
in this system. We begin with the Hubbard Hamiltonian, 



H — Hkin + Hi n t ) 



with, 



and, 
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The Hamiltonian contains a kinetic part Hkm with a hop- 
ping integral from site j to i, and an interaction part 
with a local Coulomb repulsion U for electrons on the 
same site. c\ a (cj a ) is the creation (annihilation) oper- 
ator of an electron at site i with up or down spin a. 
rii a = c\ a Ci a measures the number (0 or 1) of electron 
at site i with spin a. The Hamiltonian, Eq. [T] con- 
tains the key ingredients for correlated up and down spin 
electrons on a lattice: the competition between delocal- 
ization of electrons by hoping and their localization by 
the interaction. It is one of the mostly used models to 
study the electronic correlations in solids. 

In the absence of the interaction U, the ground state is 
characterized by the Slater determinant comprising the 
Hartree-like wave functions (HWF) of the uncorrelated 
electrons, \4>o)- When U is turned on, the weight of the 
doubly occupied sites will be reduced because they cost 



2 



an additional energy U per site. Accordingly, the trial 
Gutzwiller wave function (GWF) \ipa~) is built from the 
HWF |Vo), 

|^ G )=g D |Vo>. (4) 

The role of g D is to reduce the weight of the configura- 
tions with doubly occupied sites, where D = J2i n i1 n il 
measures the number of double occupations and g(< 1) 
is a variational parameter. In fact, this method corrects 
the mean field (Hartree) approach for which up and down 
spin electrons are independent, and, some how, overesti- 
mates configurations with double occupied sites. Using 
the Rayleigh-Ritz principle, this parameter is determined 
by minimization of the energy in the Gutzwiller state 
\ipG ) , giving an upper bound to the true unknown ground 
state energy of H . Note that to enable this calculation 
to be tractable, it is necessary to use the Gutzwiller's 
approximation which assumes that all configurations in 
the HWF have the same weight. 

Nozieres [14j | proposed an alternative way which shows 
that the Gutzwiller approach is equivalent to the renor- 
malization of the density matrix in the GWF. It can be 
formalized as 

PG = TVoT. (5) 
The density matrices pg = \iPg)('4'g\ and po = \ipo)(4>o\ 



are projectors on the GWF and HWF respectively. T is 
an operator which is diagonal in the configuration basis; 
T = UiTi where Tj is a diagonal operator acting on site 
i, 



T i \L i ,L') = J^-fL\L i ,L>). (6) 

Here, Li is an atomic configuration of the site i, with 
probability p(Li) in the GWF and po(Li) in the HWF re- 
spectively, whereas L' is a configuration of the remaining 
sites of the lattice. Note that this prescription does not 
change the phase of the wave function as the eigenvalues 
of the operators Tj are real. The correlations are local, 
and the configuration probabilities for different sites are 
independent. 

The expectation value of the Hamiltonian is given by, 

(H) G = Tr(p G H). (7) 

The mean value of the on-site operators is exactly cal- 
culated with the double occupancy probability, di = 
{n^Tii^G- di is the new variational parameters replacing 
g. Using Eqs. (|5I6[) , the two-sites operator contribution 
of the kinetic energy can be written as, 



Po(L' a ) 



(8) 



where L' a and L a are the only two configurations of spin 
a at sties i and j that give non-zero matrix element to the 
operator in the brackets. The summation is performed 
over the configurations of opposite spin L_ a . The prob- 
abilities po in the HWF depend only on the number of 
electrons, whereas the p in the GWF also depends on di. 

After some elementary algebra, one can show that the 
Gutzwiller mean value can be factorized into, 

( c L c j<?)g = Vq^(Ci a c j<T )oy/q—, (9) 

where these renormalization factors q^ are local and can 
be expressed as, 

- riia - ni- a + di + \fd i )^Jn l _ a - d l 
\Jn l(J {\ - n i(J ) 

(10) 

We used shorthand notation for (TJj CT ), the aver- 

age number of electrons on the considered "orbital-spin" 
in the HWF, which could be site and/or spin indepen- 
dent if the state is homogeneous and/or paramagnetic; 
this is the case we consider here. 



As seen in Eq. [SJ the kinetic energy of the non- 
interacting electrons state, (c| .c J - tr )o = £ kin> * s renormal- 
ized by a factor of q which is smaller than one in the cor- 
related state, and equal to one in the HWF. This factor 
can be interpreted as a direct measure of the correlation 
effect. Indeed Vollhardt [l5| has shown that l/q= m* jra 
where m* is the effective mass and m is the bare mass 
of the electron. Thus a q close to 1 corresponds to less 
correlated electron system and smaller q value reflects 
enhancement of the correlation effect. Eq. [7| leads to the 
variational energy, 

E(d) = (H) G =2qe° kin + Ud, (11) 

which can be minimized with respect to the variational 
parameter d. In the case of half filling (n=l/2), we can 
show that the minimization condition is given by, 

'-^-isb- ,12) 
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and, 



q=(l- 



U 



(i64 m ) ; 



•)• 



(13) 



It implies that if the Coulomb repulsion U exceeds a 
critical value U c = 16e° in , q is equal to zero, leading 
to an infinite quasiparticle mass with a Mott-Hubbard 
Metal-Insulator transition which is also known as "the 
Brinkmann-Rice transition" [16| . as these authors first 
applied the Gutzwiller approximation to the Metal- 
Insulator transition. 

To determine the magnitude of U c we need to calculate 
the kinetic energy of the p z -electrons in graphene. We use 
the recursion method to first calculate the local density 
of states (LDOS), i.e., projected on a given orbital: 



N^{E) =Y,KE ~ E n ) |<Vi I 



ir ij^o h + ir\ — ti 



(14) 



where \n) are the eigenstates associated with eigenvalues 
E n , and using it we calculate the kinetic energy. Let us 
briefly remind the recursion method of finding the density 
of states. We follow Ref. [l7| in describing the recursion 
method. 

For a given normalized state in which we want to 
calculate the LDOS, we can always associate a recursion 
basis which is constructed by a Schmidt orthogonaliza- 
tion procedure starting from the set of states |V'i)i H \ipi), 

H 2 \^i),...,H n ~ 1 \i{ji) Let us consider H\ip-\). We can 

decompose it in a component parallel to \tpi) and a com- 
ponent orthogonal to So we can write, 



H\ip 1 ) = a 1 \i; 1 ) + b 1 \Tp 2 ). 



(15) 



If the coefficient b\ is chosen real and positive then \tp2) 
can be defined in a unique way. Consider then H\ip2) , this 
vector can be decomposed in a component parallel to the 
space spanned by \ip\), \ip2) and a component orthogonal 
to this space. Then we get, 



H\$ 2 ) = a 2 |V>2) + ft'il^i} + hlfa 



(16) 



Furthermore, since H is a hermitian operator we deduce 
that b[ — b\. If the coefficient bi is chosen real and 
positive then ^3) can also be defined in a unique way. 
We can repeat the process leading to the construction of 
a set of states \ip n ) which are orthonormal and satisfy, 



H\tp n ) = a n \ip n ) + b n -i\ip n -i) + b n \ip n+ i). 



(17) 



In the basis of \ip n ) the Hamiltonian is tridiagonal. An 
important property of the sates \ip n ) is that they spread 
progressively from the initial state. For example if \ipi) 
is located on an atomic orbital p z then the state \tp n ) can 
have non-zero components only on orbital that is reached 
by n — 1 applications of the Hamiltonian H, i.e., n — 1 
hopping from initial orbital. 



KE and DOS as a fucntion of filling 
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FIG. 1: (Color online) The density of states (DOS) and the 
kinetic energy (KE) as a function of filling factor in graphene. 
Energy is measured in unit of t. 



The recursion basis establishes an order between states 
that are attained at shorter or longer times when the ini- 
tial states is Thus, roughly speaking one expects 
that it is important to treat properly the states with 
smallest indices n, and in general it will have smaller 
effect if states with larger indices are approximated or 
simply not considered. 

From a technical point of view the interest of the recur- 
sion basis is that it can be used to calculate vectors like 
f(H)\i[)i) where f(H) is any operator that is a function 
of the Hamiltonian H . This vector can be decomposed 
in the basis of the states \ip n ) as, 



/(ff)|Vi>=l>»><iM/(Wi 



(18) 



When one applies this method to a function, f(H) = 
(z — one can show that, 



R(z) = (ipi\ 



z- H 



1 



z — ai — 



( 19 ) 

In the simplest case where the spectrum does not have 
a gap, the coefficients a n and b n tend to a limit and 
boo respectively. If these limits are obtained to a good 
precision one can write 



R(z) 



where, 



(20) 



z — a\ 



- b n R T( z ) 



Rt(z) 



b 2 



(21) 
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FIG. 2: The kinetic energy renormalization factor as a func- 
tion of on-site Coulomb potential. For practical values of po- 
tential the renormalization of the kinetic energy is very small. 



is known as the terminator of the series. This infinite 
expansion means that Rt{z) obeys a simple relation, 



Rt(z) 



from which one deduces, 



1 



Rt(z) = 



-a- b 2 R T (z) ' 



±y/(z-a) 2 -4b 2 
262 ' 



(22) 



(23) 



So Rt(z) is the resolvent for a semi-elliptic density of 
states with energies in the range a — 2b < E < a + 2b. 
In the case of graphene, due to the the symmetry of the 
honeycomb lattice, all a n = 0. On another hand the 
presence of the pseudo gap in the middle of the band 
leads to two alternating limits for b dd and b even . Thus 
we use a second order terminator, which can be written 
as. 



Rt(z) = 



b 2 



(24) 



1 Rt(z) 



which can be easily solved. 

This continued fraction expansion of the diagonal el- 
ement of the resolvent provides an efficient way of com- 
puting the LDOS using the relation, 



N{E) = — -ImR(E). 

7T 



Once we get the density of states as a function of filling 
factor, we can calculate the kinetic energy per spin using, 



EN(E)dE. 



(26) 



a-2b 



(25) 



The calculated density of states and kinetic energy as 
a function filling factor is presented in Fig. 1. The en- 
ergy is measured in terms of t = 3.0eV(— ppw), the 
nearest neighbor hopping energy. From the figure we 
can see that the kinetic energy per spin at half filling, 
£ kin — 2.4eU. So the critical on site repulsion required 
for the Brinkmann-Rice transition is Uc = 38.4eU which 
is far from being a plausible magnitude of the Coulomb 
interaction in graphene. 

We also calculate the kinetic energy renormalization 
factor q as a function of the on site Coulomb energy. 
The result is shown in Fig. 2. We can clearly see form 
the figure that the renormalization of the kinetic energy 
due to the on-site Coulomb interaction is very small for 
practical values of the interaction. For example, even 
for the interaction as high as lOeV the kinetic energy is 
reduced by only about 5%. 

In conclusion, we used Gutzwillcr method to study the 
Brinkmann-Rice transition in single layer graphene. We 
calculate the density of states and the kinetic energy 
of electrons at half filling using recursion method. We 
find that the kinetic energy of electrons in graphene is 
very high such that the potential energy required to have 
the metal insulator transition is unpractically high. The 
on-site Coulomb interaction has very small effect on the 
renormalization of the kinetic energy. 
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